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It is shown that the limit t — t' — » oo of the equilibrium dynamic self-energy can be computed 
from the n — ► 1 limit of the static self-energy of a n-times replicated system with one step replica 
symmetry breaking structure. It is also shown that the Dyson equation of the replicated system 
leads in the n — » 1 limit to the bifurcation equation for the glass ergodicity breaking parameter 
computed from dynamics. The equivalence of the replica formalism to the long time limit of the 
equilibrium relaxation dynamics is proved to all orders in perturbation for a scalar theory. 
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I. INTRODUCTION 

Spin glass and structural glasses are characterized by the presence of a complex structure of stable and metastable 
states p}. The basic idea @, y, 0, la, la, 0, @] , supported by the behavior of some mean- field spin glass models, is that 
the glassy behavior arises because of the emergence of an exponentially large number of metastable states that breaks 
the ergodicity and prevents the system from reaching the true thermodynamic equilibrium state. The logarithm of 
the number of metastable state is commonly called "complexity" or "configurational" entropy. The glass transition is 
then associated with the emergence of a non-zero configurational entropy. By extending the replica method, originally 
developed for disordered systems, the occurrence of a non-zero configurational entropy can be conveniently studied 
within the replica formalism Q. 

The basic idea of the replica approach 0, [ly 1 , [TJ, [l2|, [H, H3] is to consider the equilibrium a thermodynamics of 
m copies of the original system interacting among them via an infinitesimal attractive coupling. If the free-energy 
landscape breaks down into (exponentially) many metastable states, the m copies will then condensate into the same 
metastable state since their relative distances in such states are smaller than in the liquid (paramagnetic) phase. The 
replica method allow for an equilibrium analysis of the properties of the metastable states, in spite of the fact that 
these are originally defined in a dynamic framework. In particular the properties of the original system is recovered 
continuing the replica number m to m < 1 and taking in the limit m — * 1. 

If one applies the present method to mean-field p-spin spin glass models, whose low temperature phase is described 
by a one replica symmetry breaking step (1RSB), the result agrees with that from a dynamical study. We observe 
that in this case the number m of replicas corresponds to the break point of the 1RSB solution of the conventional 
replica approach, and hence the m — * 1 limit clearly corresponds to the onset of the glass phase. 

All checks of the equivalence between the replica approach and the dynamical approach we are aware of 15, 1G, 
E3, Il8j . use some approximations to compute the relevant correlation functions of the many-body problem. These 
leave open the possibility that the equivalence is just a consequence of the limited number of dynamical diagrams 
considered in the various approximations. 

To our knowledge the equivalence of two approaches for a general system is not yet proved. The proof of this 
statement is the main result reported in this paper. To be more specific, in this paper we show under rather general 
assumptions, that the long time limit t — t' — > oo of the equilibrium two-times correlation function of a glassy system 
can be computed from the n — ► 1 limit of the static theory of an n-times replicated system with a 1RSB structure. 
The proof is based, following similar works |l9l . |2C| on the connection between static and equilibrium dynamics, on a 
suitable summation of classes of dynamical diagrams of the dynamic perturbation theory generated by the Langevin 
dynamics, and showing that in the limit t — t' — > oo these approaches those of the n — ► 1 limit of the equilibrium 
theory of an n-times replicated system with a 1RSB structure. 

The reader will not find in this paper the physical consequences that can be obtained for specific s yste ms or models 
by using the results reported here. For these the interested reader is referred to, e.g., Ref. flol IT^.TT7l Il8j or to the 
forthcoming paper (2lT | where these will be applied to a simple toy system. 



'Electronic address: andrea.crisanti@romal.infn.it 



2 



To keep the notation as simple as possible we shall consider only scalar fields. The generalization to more complex 
fields, e.g., vector field, is straightforward. Our starting point is hence the Langevin equation of the form 

dtp(t) 5H[<p] 

which describes the purely relaxation dynamics towards equilibrium of the scalar stochastic field tp in presence of the 
stochastic force rj. In general both tp and rj can be functions of both time and space coordinates however, since we 
are interested into the time behavior of correlations, space coordinates can be safely neglected. As a consequence we 
shall drop any explicit space dependence of fields. The Hamiltonian H[ip] governs the behavior of the system. Here, 
when needed to illustrate the calculations, we shall use 

HM = r -v 2 + jy (2) 

which describes a zero-dimensional scalar tp 4 theory. 

The stochastic force r\ has a Gaussian distribution of zero mean and second moment fixed by the Einstein's relation 

( V (t)v(t')) =2TS(t-t') (3) 

where T is the temperature. Indeed, as can be verified by means of the associated Fokker-Planck equation [22j], eq. 
© guarantees that the time dependent probability distribution function P[tp, t] eventually converge for t — > oo to the 
equilibrium distribution 

P cq [tp] cxexp(-/3ffM), (4) 

where (3 = T 
Given an i 

generated by the Langevin equation (fl]) can be written as the path integral [23j, |24j, |25|, |26|, |27|, [2 

(O) = [ VtpVtp 0[ip}e- s[v ^ ] (5) 



Given an initial condition tp(to) the expectation value of a generic observable Q\tp] over the stochastic process 

IMF 



where 

SH[p] 



<p]= dt < tp 



la 



d t tp 



tp> (6) 



The dynamical functional S[<p, tp] takes the form of a statistical field theory with two independent sets of fields, the 
original field tp and the response field t/3, so that all well established machinery of statistical and quantum field theory 
can be applied. The field ip is called "response field" since it is conjugated to an external field h. As a consequence 
arbitrary response functions can be generated by taking correlators involving ^-fields. In particular the dynamic 
susceptibility or response function G(t, s) reads 



oft.! -'<*'» 



Sh(s) 



= t>s. (7) 

h=0 



From now on we absorb the factor (3 into the definition of the response function. 

We note that in deriving the dynamical functional S[ip, tp] we have assumed the Ito prescription for the stochastic 
calculus. This implies that G(t, t) = 0. Accordingly, in a perturbative expansion of averages ([5]) all diagrams with at 
least one loop of the response function can be neglected. 

In equilibrium the response function G(t, s) is related to the two point correlation function C(t, s) = (tp(t) tp(s)) by 
the fluctuation-dissipation theorem (FDT): 

G(t, s) = 9(t - s)— C(t, s) - -6{t - s)— C(t, s) (8) 

We note that time translation invariance of equilibrium implies that all two-points correlators are function of the time 
difference only. This functional dependence is always assumed, even if not explicitly shown, every times we consider 
equilibrium correlators. 
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In equilibrium the calculation of the two point correlation function C(t, t') can be reduced to that of the dynamic 
self-energy E<^(i,t'). Indeed by using standard methods of statistical field theory [2^| and the FDT relation © one 
obtains for C(t,t') the following equation, valid for t > t': 



[d t + H"{0] - t)] C(t, t') + ds E^(t, s) 8 s C(s, t')=0 t>t' 



(9) 



where 



H"[0] 



PH[<p] 



Sip 2 



(10) 



is the quadratic part of the Hamiltonian H[tp\. For simplicity we assume that in equilibrium (ip) = 0. 
In glassy systems the correlation function C(t,t') does not vanish in long-time limit t — t' — > oo: 

lim C(t,t') = C(oo) ^ 



(11) 



signaling the breaking of the ergodicity [30]. The value of the "ergodicity breaking parameter" C(oo) is obtained by 
taking t — t' — > oo limit of eq. ([9]) : 



[H"[0] - E^(0)] C(oo) + [C(oo) - C(0)] E^(c 







(12) 



where E^(oo) is the t — t 1 — > oo limit of E^^(t, t'), while C(0) = C(t,t) and E^(0) = E^^(t, t) are the equal-time 
value of the correlation and self-energy. The latter can be eliminated using the relation 



[H"[0] - E^(0)] C(0) = 1 

that follows form the t — > t' + limit of equation ([5]) and lim t _ ><: <+ dtC(t,t') 
cq. (|13p we finally obtain the bifurcation equation 



C(oo) 



C(0) E^(co) 



C(oo) 



C(0) 



(13) 

T. By combining together eq. (fl^)) and 



(14) 



In this approach the glass transition is signaled by the appearance of a non-trivial solution of the bifurcation equation. 

The calculation of the dynamic self-energy E^(i, t') is in general a non trivial task and approximations are usually 
required to deal with the diagrams of the dynamical perturbation theory. In this paper we shall show that in the limit 
of t — t' — > oo the calculation of the dynamic diagrams simplifies and the bifurcation equation (|14[) can be obtained 
from a purely static calculation of a n-times replicated system with a 1RSB structure. 

The paper is organized as follows. In Section [ill we shall study the structure of the equilibrium diagrammatic 
expansion of the dynamic self-energy T,^(t,t'). We shall show that for t > t' E^(t, t') can be decomposed into the 
sum of classes of diagrams obtained by a suitable rearrangement and partial summation of dynamical diagrams. In 
Section UTT1 we shall consider the limit t — t' — > oo of the self-energy E^(f, t'). Here, using both a sum rule approach 
and a diagrammatic approach, we shall prove that in the limit t — t' — > oo all time integrals of the diagrammatic 
expansion can be evaluated and that the self-energy S^(oo) is a function of C(0) and C(oo) only. In Section HVl we 
develop the replica calculation. We shall show that E^^(oo) evaluated in Section Hill is equal to the static self-energy 
of an n-times replicated system with a 1RSB structure when the limit n — > 1 is taken. This result allows us to derive 
the bifurcation equation (|14p on a purely static calculation. 



II. EQUILIBRIUM DYNAMIC SELF-ENERGY: t t' FINITE 



The self-energy gives the correction to the free theory correlators when interactions are taken into account. Its 
calculation is in general rather difficult and one is obliged to use some approximation schemes. Following standard 
field theory procedures [1^, [3l| a systematic perturbative calculation of the self-energy can be established in terms of 
the so called proper vertex functions. Diagrammatically these quantities are represented by the sets of one-particle 
irreducible (1PI) diagrams, i.e., by diagrams that do not split into two subdiagrams by cutting any single propagator 
line, with all the external incoming ip and outgoing ip legs amputated. 
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A. Equilibrium Dynamic Diagrams 

The self-energy E^ (£,£') is related to the proper vertex with the two external outgoing ip legs removed. The 
equilibrium diagrammatic expansion of the dynamic self-energy E^(t, i') takes then the following form 



s^(m') = -CfjD— + -CQD^ + - as) 

where each line connecting the two vertices represents the full correlation function [32|, [H, [34[ 



C( S , S ') = {<p(s)<p(s')) = (16) 

s s 

The full (left) vertex with 1 + r external legs is the sum of all connected 1PI dynamic diagrams with the outgoing ip 
leg at i, the wiggly line in (|15[) . and r incoming ip legs at times si, . . . , s r (< t) removed. To each internal y^-line 
connecting the internal vertices at times s' and s is associated the full correlation function C(s, s') while to each 
internal tptp-line connecting the internal vertices at times s' and s is associated the full response function (32l . l33l . l34l ] 



G(s, s') = {tp(s) tp(s')} = 'wwv (17) 

s s' 

with s > s'. All internal times are integrated from t' to t. 

The empty (right) vertex with 1 + r external legs is built from the same 1PI diagrams of the full vertex but with 
all external times times t, s± . . . , s r equal to t! . As a consequence, since in equilibrium the correlation and response 
functions are related by FDT, it is given by the topological equivalent vertex obtained from the associated static 
equilibrium theory described by the canonical distribution ^ [2fJ. 

The structure just described can be understood as follows. First of all we note that since closed loops of response 
function vanish, each internal vertex of the 1PI diagrams contributing to the self-energy E^(i, i') is connected via 
</?<£>-lines to one or the other of the external wiggly lines, but not to both. This means that the diagram can be divided 
into two subdiagrams made of all the vertices connected to the same external wiggly line. The two subdiagrams are 
clearly 1PI and joined together by 2 or more ipip-lmes. It is easy to recognize the two subdiagrams as the full/empty 
vertices. 

Alternatively one can invoke the general diagrammatic expansion of the correlation function (ip(t) <p(t')) as the 
matching of two tree expansions [35j |. one for ip{t) and one for f(t'), and note that the lines connecting the full/empty 
vertices are the lines joining the two tree expansions. We note that the assumption that the system is at equilibrium 
implies that all two-times quantities depend only on time difference. Then the contribution from the tree expansion 
of <p(t') averaged over noise and equilibrium initial conditions, the empty vertices, cannot depend on time t' and must 
be equal to that obtained from the equilibrium static theory described by the canonical distribution (0| with all lines 
equal to the static equilibrium correlation function C(t! , t 1 ) = C(0). 

The internal structure of the empty vertices can be inferred by using the following dynamical functional, see 
Appendix [Cl 

/>oo 

S'[<p,<p] = S[<p,<p]+ dsH[<p( s )]6{s-t ) (18) 

to impose statistical equilibrium with the canonical distribution Q at the initial time to . The analysis of the dynamic 
diagrams generated by S'[ip, tp] shows that the effect of the last term is that of canceling out from all dynamic diagrams 
for E^(i, t') the contribution from times to < s < t' yielding for E^^i') the equilibrium dynamic diagrammatic 
expansion discussed above. This can be understood on a general ground as follows. The canonical distribution ^ is 
a stationary solution of the associated Fokker-Planck equation, and hence the probability distribution of tp remains 
canonical for any time past to. Consequence of this is that all quantities evaluated from S'[(p, ip] cannot depend on to- 
This guarantees, for example, that all two-times quantities depend only on time differences. In evaluating T,^(t,t') 
for t > t' we can then choose for to in (fT5|) any value < t' . The invariance property ensures that we always get the 
same diagrams. Clearly the simplest choice is to = t' which, in turn, implies that in the diagrammatic dynamical 
perturbative expansion all free times must be integrated from t' . We have seen that the dynamical diagrams for 
E^(i,t') generated by S[ip, tp] can be divided into two subdiagrams, joined by r > 2 correlation lines, by grouping 
together all vertices connected by response c^-lines to the same external (amputated) wiggly line at t or t', respectively. 
The dynamical diagrams generated by S'[tp, ip] can be divided in a similar way into two subdiagrams connected by 
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s 




t t' 
(b) 



t t' 
(c) 



FIG. 1: Third order dynamic diagrams with equilibrium condition at to. 



r correlation lines just grouping together all vertices connected by response tptp-lmes to the external (amputated) 
wiggly line at t. It is easy to realize that this procedure leads to the same full vertex obtained from S[tp, ip], while the 
(putative) empty vertex, the one connected to external (amputated) leg at t', contains now only contributions from 
the last term in (|18p since the equal time response function vanishes. Stated in a different way, the empty vertex 
contains only the diagrams generated by the equilibrium initial condition, i.e., it is equal to topological equivalent 1PI 
diagram with 1 + r external (amputated) legs of the associated static theory described by the canonical distribution 
(HI with all ipip-\mes equal to the static equilibrium correlation function C(t,t) = C(0). The same conclusion can be 
obtained by first dividing the dynamical diagrams generated by S'[ip, <p\ as described above, and then taking the limit 
to — > t' directly on diagrams. 

The definition of empty vertex given above uses S[<p, <p] and follows from the observation that the equilibrium 
FDT relation between response and correlation function guarantees that setting all external times of a dynamic 
diagram equal to each other reduces the dynamic diagram to the topological equivalent diagram of the associated 
static equilibrium theory [20l ]. 

We can then summarize the rules for writing down the equilibrium dynamic diagrams for Yi^{t,t'): 

1. Write down the dynamic diagrams generated by the dynamical functional S[cp, tp] using the standard dynamical 
rules, neglecting all numerical symmetry factors. 

2. In each diagram remove the minimal number of G-lines needed to divide the diagram into two disjoint sub- 
diagrams so that in the first (left) sub-diagram all vertices are connected to time t through G-lines, while in the 
second (right) sub-diagram are connected through G-lines to time t'. 

3. In the second (right) sub-diagram, the one connected to t 1 , replace all response G-lines by correlation G-lines, 
and set all time variables to t'. 

4. If after replacement two or more different dynamic diagrams lead to the same diagram count the latter only 
once. 

5. Multiply each diagram so obtained by the appropriate numerical symmetry factor and evaluate it with the usual 
rules integrating all left internal times from t' to t. 

To illustrate the above rules consider the third order dynamic diagrams shown in Fig. [1] generated by the dynamical 
functional S'[<p,<p] ([18f for the scalar zero-dimensional ip 4 theory ([2]). 

By using the standard dynamic rules the contribution of these diagrams is 



2$(t,lf) = \ f d s G(t,s)C{t,s)C 2 (t',s)C{t,t') + \ f dsC 2 (t,s)G(t', s )C(t',s)C(t,t') 




+-C 2 (tM)C 2 {t'M)C{t,t l ) 



(19) 



With the help of the FDT relation {HJ the contribution can be rewritten as 





(20) 
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FIG. 2: Third order equilibrium dynamic diagrams. 




FIG. 3: Third order base diagram. 



The relevant third order diagrams for the equilibrium dynamic rules, point 1), are the first two diagram in Fig. [T] 
The decomposition of points 2)-4) leads to the equilibrium diagrams shown in Fig. [5] The two diagrams are easily 
evaluated and one recovers the result ([20)1 . 



By reversing the rules to write down the equilibrium dynamic diagrams it follows that each equilibrium dynamic 
diagram can be obtained by a suitable decoration of a base diagram, i.e., of the diagram with the same topology of 
the equilibrium dynamic diagram. For example the base diagram for the dynamic diagrams of the previous example 
is the one shown in Fig. [3l 

The equilibrium dynamic diagram is generated from the base diagram by: 

1. dividing the base diagram into two sub-diagrams by cutting some lines to reproduce the topology of point 2) of 
the equilibrium dynamic rules. 

2. In the (left) sub-diagram associated to the full vertex attach a 0-line to each vertex to generate the desired 
G-line connection structure. 

It is clear that the base diagram is by construction the dynamic diagram in which all G-lines are replaced by C-lines. 
As a consequence the self-energy base diagrams are the static self-energy diagrams of the associated equilibrium static 
theory described by the canonical distribution ((3]). 



From the structure (|15p of the diagrammatic expansion of the self-energy it is easy to realize that also each dynamic 
diagram contributing to the full vertex with 1 + r external legs can be obtained from a suitable (static) base diagram 
by attaching to each internal vertices one 0-line to generate the desired G-line structure. Figure 2] shows a base 
diagram for the 1 + 3 full vertex and the three different dynamic diagrams that can be generated. 

A convenient way of imposing that to each vertex is attached one and only one G-line is by means of anti-commuting 
Grassman variables. Following Ref. |20l | we then introduce for each vertex the pair of conjugated Grassman variables 
(a a ,a s ), where s is the time variable label of the vertex. By using the FDT relation {8J and the properties of Grassman 
variables it is easy to sec that 



B. Self-Energy Base Diagrams 



C. Partial summation of full vertex diagrams 



C[s — s' + a s a s 9(s' — s) — a s >a s i0(s — s')] 



C(s — s) + a s a s 6(s — s)d s > - s C(s — s) + a S 'a s >6{s — s')9 s <_ s G(s — s) 
C(s', s) + a s a s G(s', s) + a s >a s >G(s, s'). 



(21) 
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The last equality follows from even parity of the correlation function C(t, s) = C(s, t). Equation (|2Tj) has the following 
simple diagrammatic representation: 



C[s — s + a s a s 9(s' — s) — a s ia S '9(s — s')] = + a s a s -ww + a s ra s : wwv (22) 

s s' s s' s s' 

Consider now a 1PI base diagram B for the 1 + r full vertex made of J internal bare vertices, i.e., vertices without 

external legs, and 1 + E external bare vertices, i.e., vertices with at least one external leg: 



B := 1 — ( I V- E (23) 




Each bare external vertex has m = 1,2, . . . external legs, with X^fLo Tit = 1 + r. Then from (|22p and the properties 

dada aa = 1, da 1 = da 1 = (24) 



of Grassman variables it follows that the sum Ab(sq, s r ) of all dynamic diagrams of the 1 + r full vertex that can 
be generated from the base diagram B can be obtained as: 

1. Order the 1 + E + I bare vertices of B so that: the label corresponds to the external vertex attached to the 
out-going (p-leg at time t, the labels 1, .... i? to the remaining external vertices and labels E + 1, . . . , E + 1 to 
the internal vertices. 

2. Assign to vertex the time variable so = t, while to each vertex labeled by i = 1, . . . , E + I assign the time 
variable Sj and the pair of conjugate Grassman variables (a Si ,a Si ). From the time-oriented nature of the 
dynamic diagram it follows that sq > Si > t' for i = 1, . . . , E + I. 

3. Assign to each line connecting the vertices i and j, with i,j ^ 0, the function 

C[si - Sj +a Si a Si 9(sj - Si) -a Sj a Sj 6{si - Sj)] (25) 
and to each line connected to the vertex the function 

C[si - s + a Si a Si 6(s - s,;)] (26) 

with i,j = 1,...,E + I. 

4. Integrate over all Grassman variables, 



;.=i 

to ensure that all vertices have one, and only one, <^-line attached to them 



. E+I 

J [] da Si da Si ••■ (27) 




wu^^^- C(t,t') -^^^^ww 

t t' 



FIG. 5: Factorization of generic dynamic self-energy diagram for t — t' — » oo. 



5. Integrate all internal time variables s$ from i' to sq = t, 

rsa E+I 



r-s -^TJ 

j n d ^ 



i=E+l 

When these steps are translated into formulae we end up with: 

„ SQ E+I „ E+I E+I 

A B (so,s 1: . . . ,s E ) = M(B) I Yl ds i / H da Sz da Si ] [ K 0i C[s l - s 

i=B+l J i=l i=l 

x II II ^i*^ 3 ' _ S J +a s% a s% 9(s :j - -a Sj a Sj 9(si - Sj)] n ** 

3 = 1 i=l 



a Si a Si U(s - Si 



/sq E+I E+I 
]J <ls, / ] [ da Sz da Sl 
i=E+l J »=1 



»=J3+1 

x IT II ~ S J + ^siCts^isj - - a Sj a Si 9(si - sj 

j=0 »=i 



where M(B) is symmetry factor of the base diagram £? and 

1 if vertices i, j are directly connected 



otherwise 



(28) 



(29) 



riij = riji — # of lines connecting vertices i and j 

give the connection topology of B. The last equality in eq. (|29]) follows from the observation that #(s, — sq) = for 
all i = 1, ...,£' + /. 

It is easy to verify that the integration over Grassman variables produces all possible dynamic diagrams, with the 
correct weighting factor, that can be generated from the base diagram B. 



III. EQUILIBRIUM DYNAMIC SELF-ENERGY: THE LIMIT t - t' -> oo 

In the limit t — t' — > oo some simplihcations occur in the calculation of the dynamic self-energy diagrams. Each 
bare vertex making up the 1 + r full vertex is connected to the bare vertex attached to the external <^-leg at time t 
through G-lines then, since 



G(s, s') = 9(s - s') d s ,C{s, s') ^ iff s' ~ s, 



(30) 



it follows that the time variables s%, . . . , s r of the remaining r external legs are Sj ~ t ^> t' for t — t' — > oo. As a 
consequence in the diagrammatic expansion of the self-energy E^(f , t 1 ) we can replace for t — t' — > oo all correlation 
lines connecting the full and empty vertex pair with C(t,t') = C(oo) and the generic diagram factorizes as shown in 
Fig. [5j The overbar on the 1 + r full vertex means that all external time variables si, . . . , s r are integrated from t' to 
t. 

The integration in the 1 + r full vertex can be done considering first the contribution of all dynamic diagrams that 
can be generated from a 1 + r full vertex base diagram B and then summing up the contributions from all possible 
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base diagrams. Neglecting the full-empty vertex connection symmetry factor, the contribution of all dynamic diagram 
generated by the 1 + r full vertex base diagram B with 1 + E external and / internal vertices, can be written as: 

f J] d Si A B (t, Sl , . . . , s E ) n C( Si , t') ni C(t, t'T - 1 A (0) (f ) (31) 

Jt ' i=l i=l 

where rij is the number of external legs of the z-th external bare vertex of B, and A^°\t') = A'°)(C(0)) is the 
contribution of the 1 + r empty vertex. [36[ 

In the limit t — t' — > oo the first term of eq. (|3"Tj) is not zero only if all Sj ~ t 3> t' and we can replace all C(si, t') 
in the second term with C(t,t') = C(oo). Equation (|3"Tj) then factorizes, cfr. Fig. as 

A B (i)C(i,i7A (0) (i') (32) 

where 

As(i) = / TTdsi A B (s , si, . . . , Sfi)| So=t ■ (33) 
i=i 

Finally by adding the appropriate full-empty vertex connection symmetry factor and summing eq. (|32p over all 
possible 1 + r base diagrams B one recovers the complete dynamic diagram contribution in the limit t — t' — > oo from 
the 1 + 7- vertex to X^(<, t'). 

A. 1 + r full vertex base diagram 

From the definition of As(so) (1331) and the expression (f2"9")) it follows that 

So JV JV AT JV 

A B (s ) = M(B) JJdsi /[](iQ Si rf5 Si [][[if J1 C[s 1 - Sj +5 Si a s ,%-s 1 )-Q S] a s i(s,- Sj )p (34) 

i=l J t=l j=0i=l 

where A = / + is the total number of bare vertices making up B, excluding the final one connected to the external 
ip-leg at time s = t. This expression can be simplified by integrating over si, . . . , sn in a fixed order because from 
FDT it follows that: 

!C[x - y + a x a x 8(y - x)} if x < y 
(35) 
C[y — x + a y a y 9(x — y)\ if x > y 

Thus ordering si, . . . , spj so that: 

(si,...,s N ) -> (s h ,...,Si N ) : s tp > s tq iiq>p (36) 
the integrand of eq. ([34]) can be rewritten as 

N N 

IT IT K i l C ^ Si ~ S J + s *s a s^( s j - s*) - a Sj a Sj 6(si - Sj)] 71 ** 

j=0 i=l 

W O.N 

^ II II *<«<p C K - % +as lp a Slp 9(s tq - s ip )] n ^p (37) 

p=l q<p 

Inserting this expression into (I34[) and using the identity: 



so 



*=i (i 1 ,...,i w )eP(i,...,iV)- /t ' 
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FIG. 6: Base diagram vertex numbering. 



where P(l, . . . , N) are the N\ permutations of (1, ... , N), we end up with 



><n 

P =i 



0,N 



dSi2 



dsi 



A B (s a ) = M(B) 

(ii,...,i N )eP(i...,w) 

N r c 

da Szp da Stp J| C[s lp - s lq +a 



Us "(Si 



9<P 



(39) 



The integration over Grassman variables is now diagonal and can be performed. A straightforward algebra leads to 



A B (s ) = M{B) 



E 



ds 



(ii,...,i JV )eP(i...,ATj *' 



m(b) yi n 

(i 1 ,...,i JV )eP(i...,iv)P=i Jt 



' v 



dSi 2 



N 



d 



dsi » n q- 



P =i 



0,N 



9<P 



dsi 



p d Sl 



]\K lqlp c( Slq , Stp r^ 

.9=0 



(40) 



As simple example of eq. (HHJ) consider the base diagram B shown in Fig. [U The bare vertices of B are numbered 
as shown in Fig. [6] There are two possible orderings, namely: {012} and {021}. We use the short-hand notation 
for Sq, 1 for Si and so on. Then from eq. (T4"0")) the contribution of this base diagram is 



A B (s ) 



^J dl j d2d 1 C(0,l)d 2 [C(0,2)C(1,2) 2 ] +\f d2 J dld 2 C(0,2)d 1 [C(0, 1)C(2, 1) 



(41) 



The factor 1/2 is the symmetry factor of the base diagram B. 

In the direct calculation we have to evaluate the three dynamic diagram generated by B shown in Fig. 0] using the 
dynamic rules and FDT, and then add the results. The first diagram from left, with the vertex numbered as in Fig. 
O leads to: 



dl J d2G(0,l)G(l,2)C(l,2)C(0,2) = | / dl / d2d x C{Q, 1)9 2 C(1, 2) 2 C(0, 2) 



while the second to 

r0 



d2j^ dl G(0, 2)G{2, 1)C(2, 1)C(0, 1) = - J d2 J dl <9 2 G(0, 2)d 1 C(2, 1) 2 C(0, 1) 



and finally the third to 

r0 r 



-J dl J d2G(0,l)G(0,2)G(2,l) 2 = -J dl J d2d 1 C{0,l 



)5 2 C(0 I 2)C(2,1) 



(42) 



(43) 



(44) 



The factor 1/2 in the third diagram follows from the symmetry of the diagram. By adding the three contributions, 
and using the identity 



M M pi pO r-2 

dl d2= dl d2+ d2 dl 
t> Jt' Jt' Jt' Jt' Jt' 



(45) 
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one easily recovers the result (1411) . 

We note that all integrations in eq. (|4"Tj) can be done. Indeed the two integrals in eq. (|4"Tj) can be combined together 
to give 



A B (s ) = J dl J d2 9 X <7(0, 1) 8 2 [<7(0, 2)<7(1, 2) : 



and performing the integral over 2 and then over 1 one ends up with 

A B (s ) = i [C(0) 4 - C*(0) 2 C(oo) 2 - 2 (7(0) C(oo) 3 + 2 C(oo) 4 



(46) 



(47) 



where (7(0) = C(t,t) and <7(oc) = C(t,t'). 

The possibility of carrying out all integrals over Sj is not a properties of this special example, but it is a general 
result valid for any base diagram B, as shown in the next subsections. 



B. Recursive integration: sum rules approach 
Equation (|4T))) can be written as 

I S lN } 



(n,.. .,%)eP(i,...,iV) 



(48) 



where 



Fb{sq\ Si i: 



p=l 



i = n /. d Sip — 



p-i 



,9=0 



is contribution of the base diagram B when the time variables s\, . . . , Sn are ordered as (s^ . 



(49) 



see eq. (|56"]l. 



We shall denote this ordered base diagram by B{sq, s^ , . 



•} 



The sum (J48J) can be performed by first selecting all terms where the pair (in-i, in) is equal to (h, k) or (k, h), and 
then summing over all possible distinct index pairs h, k. If we do this we end up with 



A B (s ) = J2 

(h,k) (i 1 ,...,i N _ 2 )eP[(l...,iV)-(h,fc)] 



E 



Fb{sq; Sij, . 



, Sh, s k } + F B {s ; s ii: . 



, Sfc, Sft} 



(50) 



where the first sum is over all distinct pairs (h, k) with h, k = 1, . . . , N, while the second sun over is all permutations 
of (1, . . . , N) once h and k are taken out. 

From eq. (|49[) it is easy to realize that the two terms in square brackets differ only for the last two integrals over 
Sh and Sk, so that we have to evaluate the sum: 



ds h d s . 



N-2 

]J K i q hC(s i<t ,S h ) 
9=0 



ds k d s 



N-2 



II K Uk C(s lq , s k ) n ^ K hk C(s h , s k ) nhk 

9=0 



+ (h ^ k) (51) 



where the second term is obtained from the first by exchanging and s k - 

By performing the s k integral in the first term and the Sh integral in the second term and rearranging the terms 
we obtain the sum rule 

F B {. ■ .,s h , s k } + F B {. ■ .,s k ,s h } = F B {. ■ -,s h \sh = s h } - F B {. . .,s h \s k = t'} - F B {. . . , s k \s h = t'} (52) 



where the dots is the short hand notation for "sq; s-ti , ■ 



of the ordered base diagram B{sq, s^, . . . , Si N _ 2 , Sh\s k — Sh} obtained from B{sq, s^, . 



The term F B {..., s^jsfc = Sh} is the contribution 



■ 5 Si 



,s h ,s k } by setting 



Sk = s h, i-e., by replacing all lines connecting vertices h and k by C(sh,Sh) — (7(0). This reduces the number 
of integrations by one since the vertices h and k "collapse" into a single composite or effective vertex. Similarly 
F B {..., Sh\s k = t'} is the contribution of the ordered base diagram B{so, s^, . . . , Si N _ 2 , Sh\s k = t'} obtained from 
B{sq, Sjj , . . . , Si N 2 , Sh, s k } by setting s k = t', i.e., by replacing all lines connecting the vertex k to any other vertex 
i 7^ k of B by C(s,,t') = <7(oo). Also in this case the number of integrations is reduced by one since the vertex k is 
"removed" from the diagram. F B {..., s k \sh — t'} is defined in a similar way. 
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This procedure can be repeated by considering triples (7, h, k) of vertices and writing 

A B (s )= E E E F B { So - Sil ,..., si,s h ,s k } (53) 

{l,h,k) (i 1 ,...,i N _ 3 )^P[(l...,N)-(l,h,k)] P(l,h,k) 

where now the first sum runs over all distinct triplets (I, h, k), the second over all possible permutations of (1, ... , N) 
once the triplet (I, h, k) has been removed and the last over the six permutations of the triplet (I, h, k). 
By using the sum rule (|52")) the last sum can be written as 

F B {...,si,s h ,s k } = 

P(l,h,k) 



(54) 



F B {.. 


■ , si, s h \s k 


= s h ] 


-F B {.. 


■ , st, s h \s k 


= t'} 


- F B {. . 


■ , Sl, Sk\s h 


=t'} 


F b {.. 


■ , Sh, Sk\si 


= s k } 


-F B {.. 


■ , Sh, Sk\si 


= t'} 


- F B {. . 


■ , s h , si\s k 


= t'} 


F b {.. 


■ , Sfc, si\s h 


= si} 


- F B {. . 


■ , Sk, S[\s h 


= t'} 


- F B {. . 


, Sk, s h \si 


=t>} 



The three term with positive sign can be combined together and give 



F B {. . . , si,s h \s k = s h } + F B {. ..,s h , s k \si = s k } + F B {. . . , s k ,si\s h = sj = 
F B {. . . , si\sh = si,Sk = s h } 
-F B {. ..,si\s h = t', s k = s h } - F B {. . . ,s h \s k = t', si = s k } - F B {. . . , s k \si = t' , s h = sj (55) 

The term F B {. . . , si\sh = si,Sk — Sh} is the contribution of the ordered diagram B{. . . , si\sh = si, Sk = Sh} obtained 
from B{. . . ,si, s h \s k — s h } by setting s h = s t , or equivalently from B{. . . , si,s h , s k } by setting s k = s h and s h = s t . 
In both cases the number of integrations is reduced by two and all lines connecting the three vertices are replaced by 
C(0). 

The term F B {. . . , si\sh = t' , s k = Sh} is the contribution of the ordered diagram B{. . . , si\sh = t' , Sk = Sh} obtained 
from B{. . . ,si, Sh\sk = Sh} by setting Sh = t' , or equivalently from B{. . . ,si, Sh, Sk} by setting Sk = Sh first and then 
Sh = t'. This means that all lines connecting the vertices h and k are replaced by C(0) while all lines connecting 
either vertex h or k with any other vertex i ^ h, k by C(oo), reducing at the same time the number of integrations 
by two. The other two terms are obtained in a similar way. 

Finally the six terms in (|54j) with negative sign can be evaluated in pairs using the sum rule (I52|) . One then gets, 
for example, 

F B {. . .,si,Sh\s k = t'} + F B {. . . , s h ,si\s k = t'} = 

F B {. . . , si\s h = si, s k = t'} - F B {. . . , si\s h = t! , Sk = t'} - F B {. . . , s h \si = t! , s k = t'} (56) 

Even if the notation should be now clear it is useful to stress that F B {. . . , si\sh = si, s k — t'} is the contribution from 
the ordered diagram B{. . . , si\sh = si, Sk = t'} obtained from B{. . . ,si, Sh\sk = t'} by setting Sh = si, or equivalently 
from B{. . . ,si, Sh, Sk} by setting Sk = t' first and only then Sh = s;. Similarly F B {. . . , si\sh — t' , Sk — if} is the 
contribution from the ordered diagram B{. . . , si \sh = t',Sk = t'} obtained from B{. . . , si, Sh\sk — t'} by setting 
Sh = t' , or equivalently from B{. . . , si, Sh, Sk} by setting Sk = t' and only then Sh = t' . The order in the construction 
is relevant since, for example, F B {. . . , si\sh — t', Sk — t'} is in general different from F B {. . . , si\sh = t' , Sk = Sh}- 
Collecting all terms we end up with the sum rule: 



P(l,h,k) 



Sl, Sh, Sk} 
























F B {... 


si\s h 


= Sl, s k 


= s h } 


















-F B {... 


si\s h 


= t', Sfc 


= s h } - 


- F B {. . 


■ , Sh\sk 


= t',Sl 


= Sk} 


-F B {. 


■ , Sk\si 


= t', s h = 


= si} 


-F B {... 


si\s h 


= Sl, s k 




F B {.. 


, si\s h -- 


= t' , s k 


=*'}+ 


F B {... 


, Sh\si = 


= t',s k = 


t'} 


-F B {... 


Sh\sk 


= s h ,si 


=*'}- 


\- F B {. . 


■ , Sh\Sk 


= t',si 




-F B {.. 


■ , s k \s h 


= t', Sl = 


t'} 


-F B {... 


Sk\si 


= Sk , Sh 


= t'}- 


VF B {.. 


■ , Sk\si 


= t', Sh 


= i'H 


-F B {.. 


■,Sl\Sk ■ 


= t! , sh = 


t'} 



(57) 

that, when inserted into eq. (|53|) . eliminates two time integrations and replaces the lines connected with the vertex 
integrated out by either C(0) or C(oo). 
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The procedure can be iterated to build sum rules for four vertices, five vertices and so. Despite the fact that the 
derivation of the sum rules for any number of vertices is straightforward we shall not push it here because these can 
be more easily obtained using the diagrammatic rules of next subsection. 

By using iteratively the sum rules and 



F B {s ; s\...} = F B {s \s = t, . . .} - F B {s \s =t',.. .} 



(58) 



all N time integrals in A b (sq) can be eliminated in turn in favor of C(0) and C(oo). This concludes the proof that 
for any base diagram B all time integrals can be performed and moreover 



A B (s ) = A B (C(0),C(oo)) 



(59) 



where A B (x, y) is a function that depends only on the topology of B. 

To illustrate the sum rule approach we conclude this subsection by reconsidering the simple example of Fig. [6l For 
this diagram N = 2, therefore from eq. (|50p and the sum rule (|52|) it follows 



A B (s ) - F B {0;1,2} + Ffl{0;2,l} 

= F B {0; 1|2 = 1} - F B {0; 1|2 = t'} - F B {0; 2|1 = if}. 



(60) 



By using now the sum rule (|58p we end up with 

A B {s Q ) = 



F B {0|1 - 0, 2 = 1} - F B {0|1 = t',2 = 1} 
F B {0|1 = 0, 2 = t'} + F B {0\1 = t',2 = t'} 
F B {0\2 = 0,1 = /} + F B {0\2 =t',l = t'}. 



(61) 



which for the diagram of Fig. [S] leads to 



A B (s ) 



C(0) 4 - C(0) 2 C(oo) 2 - 
C(0) 4 ~C(0) 2 C(oo) 2 



C(oo) 3 C(0) + C(oc) 4 - C(oo) 3 C*(0) + C(oc) 4 
-2C(0)C*(oo) 3 + 2C(oo) 4 



(62) 



One easily recognizes the result (|47|) from the direct calculation 



C. Recursive integration: diagrammatic approach 

While the integration based on sum rules described in the previous subsection can be carried on for any base diagram 
B, in practical calculations it can become quite cumbersome. Thus it would be desirable to have a simpler way of 
proceeding. In general diagrammatic methods are simpler and more transparent, for this reason in this subsection we 
present the diagrammatic approach that allow for a graphical integration directly on the base diagram B. 

The diagrammatic integration goes through the following steps: 

• Time Ordering 

The first step is to generate all possible ordered base diagrams B{so, Si lt . . . , Si N } by considering all possible 
permutations (ii, • • • , i/v) G P(l>--->Af). Once the time order of B has been fixed, each ordered diagram 
B{sq, Sjj, . . . , Si N } is decorated by orienting each line connecting two vertices by drawing an arrow in the 
direction of the time flow, i.e., pointing from the shortest to the largest time. 

The expression of F b {sq; Si 1 , . . . , Si N } is recovered by first assigning to each line of the oriented diagram 
B{sq, s^,..., Si N } the correlation function C(si p , Si q ), where i p and i q are the labels of the two vertices connected 
by the line, and then taking the product for p = 1,...,N of the derivative with respect to Sj of the product of 
the correlation functions C(sj , Sj ) associated with all outgoing oriented lines originating from vertex i p . Finally 
F b {sq; Sjj, . . . , Si N } follows by integrating the result over all Si p starting from Si N and proceeding towards in order 
[cfr. eq. (|49)) ]. In each diagram the integration starts from the vertex with the shortest time. By construction this 
vertex has only outgoing oriented lines, and there is only one of such vertex in each diagram. The integration cancels 
the arrows from all oriented lines originating from the vertex we are integrating on, and replace the original diagram 
by two new diagrams. Thus we have the following graphical integration rule. 
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BiOU) (a) (b) 

FIG. 7: Oriented base diagram B{012} and result of first diagrammatic integration. 




BI021I (c) (d) 

FIG. 8: Oriented base diagram _B{021} and result of first diagrammatic integration. 

• Time Integration 

1. In each diagram peak up the vertex with only outgoing oriented lines and replace the original diagram by 
the two diagrams obtained as follows: 

Diagram 1: Assign to the vertex with only outgoing oriented lines the next shortest time, i.e., Si p — > Si p _ 1 . 
If the vertex i p -± is not directly connected to i p draw a dotted line between the two vertices to remember 
they have the same time. Next replace all outgoing oriented lines with simple not-oriented lines drawing 
them as: full line if the line connects two vertices at different time; dashed line if the line connects two 
vertices at the same time. As before, we shall call effective vertex the group of bare vertices with equal 
time, i.e., the group of vertices connected by dashed or dotted lines. 

Diagram 2: Set the time of the vertex with only outgoing oriented lines to t' and replace all outgoing 
oriented lines of the vertex by cutted lines, i.e., replace the arrow by the "cut" sign "/". Finally 
multiply the diagram by —1 so that two new diagrams contribute with opposite signs. 

2. Merge all diagrams that differ only for the oriented outgoing lines originating from the same effective vertex 
into the single diagram obtained from any one of them by orienting all not-oriented full lines originating 
from the effective vertex according to the orientation of all diagrams we are merging to. This leaves us 
with diagrams made of only oriented lines (arrows), dashed or dotted lines and cutted lines. 

It is not difficult to realize that this rule produces diagrams where the vertex, or effective vertex, with the shortest 
time has only outgoing oriented lines. Moreover the vertices connected via cutted lines do not contribute anymore to 
the integration process. The procedure can then be iterated untill one is left with diagrams made of only dashed or 
cutted lines. The dotted lines are used for time bookkeeping and can be eliminated if not needed. 

• Value of diagrams 

The value of Ab(sq) is obtained by evaluating each final diagram with dashed lines replaced by the equal time 
correlation function C(t,t) = C(0) and cutted lines by the infinite time correlation function C(t,t') — C(oo), 
and summing up all contributions from different diagrams with the appropriate plus or minus sign. 

To illustrate the graphical integration rules we consider again the base diagram shown in Fig. [5] The two possible 
oriented diagrams £?{012} and £?{021} are shown in Figs. [7] and [8] together with the result of step 1) of the time 
integration rule. The diagrams (a) in Fig. [7] and (c) in Fig. [5] differs only for the orientation of the outgoing lines 
originating from the effective vertex made by the bare vertices 1 and 2. The two diagrams are then merged together 
into the diagram shown in Fig. [9J step 2) of the time integration rule. 
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(a) 



(0 



(e) 



FIG. 9: Merging of diagram (a) of Fig. [7] and diagram (c) of Fig. [8] 




FIG. 10: Final diagrams. 



Diagrams (b) of Fig. \7\ (d) of Fig. [5] and (e) of Fig. [5] contain only (effective) vertices with all outgoing lines or 
connected by cutted lines, thus the time integration steps 1) and 2) can be repeated for each one of these diagrams. 
This second round eliminates all oriented lines and, taking into account the —1 signs, Ab(sq) is given by the sum of 
the six diagrams shown in Fig. 1101 Evaluating these diagrams and multiplying the result by the symmetry factor 
M(B) = 1/2 of the base diagram, we and up with 



A B (so) 



C(0) 4 - C{0) 2 C{oo) 2 - C(0) C(oo) 3 + C(oo) 4 - C(0) C(oo) 3 + C(oo) 4 
C(0) 4 - C(0) 2 C(oo) - 2 C(0) C(oo) 3 + 2 C(oo) 4 
i.e., again the result (|47[) obtained previously by direct integration. 



(63) 



IV. THE REPLICA CALCULATION 



A. The 1 + r full vertex base diagram 



In the previous Section we have shown that in the limit t — t' — > oo the contribution As(t) of the 1 + r base vertex 
B to the 1 + r full vertex can be obtained by considering all possible diagrams that can be constructed from B by 
replacing the lines connecting the bare vertices by either dashed or cutted lines in all possible ways. Clearly there 
are some constraints to be considered, for example all lines connecting the same pair of vertices must be replaced 
simultaneously. In other words it not possible to have dashed and cutted line connecting the same pair of vertices, 
see e.g. Fig. [TU] 

Given any base diagram B a simple way of generating all allowed diagrams is that of attaching to each bare vertex 
i = 0, 1, . . . , N of B a label a,-. Each label can take the values 1, . . . ,n, where n is an arbitrary integer. The next step 
is that of introducing the "projector" operators S aia and (1 — 5 aia . ) and associate to each line of B connecting the 
vertex i and vertex j the "correlation" function 



C ata] = C(0)S aiaj + C(OO) (1 - 8 aiaj ) 



(64) 



By representing the projector S aiaj with a dashed line and the projector (1 — S aiaj ) with a cutted line it is easy to 
convince oneself that the sum over all Oj (i — 1, . . . , N) from 1 to n generates all possible diagrams that can be 
constructed from B. 

The interesting point of this method is that it does not only generates all possible diagrams but, as a "bonus", it 
also gives the correct weight and sign for each one of them. For a single line this is rather trivial since for n — > one 
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has the formal identity: 

/ ds — C(t, s) = C(0) - C(oo) = lim Y C ab (65) 
J t , ds n ^°t^i 

that ensures that for n — > an oriented line is not only replaced by a dashed line and a cuttcd line, but also that the 
second one comes with a negative relative sign. By using the properties of the projection operators this identity can 
be extended to a bunch of oriented lines connecting two vertices. 

In the Appendix [Al we show that the sum over labels ai with i = 1,...,N for any base diagram B reproduces in 
n — ► limit the sum rules discussed in Sec. IIIIBI This implies that Ab(sq) can be expressed in the "replica" form 

l,n N i-l 

As(t) = M(B) hm £ H H / ^ (66) 

a\ . . . .,ajv i=l J— 

We then have the following simple rule to evaluate Ab(so): 

1. Given a base diagram B assign to each bare vertex i = 0, . . . , N of B a "replica" index ai- 

2. Assign to each line connecting the vertices i and j the correlation function C ai(Li given in eq. (|64[) . 

3. Sum over all replica indexes ai with i = 1, . . . , N from 1 to n and multiply the result by the diagram symmetry 
factor M (B) . 

4. Take the limit n — > of the result. 

To illustrate the procedure, i.e., formula (|66[) . we consider again the base diagram shown in Fig. [SJ For this diagram 
eq. ([nU) reads: 

^ n n 

A~B{t) = - lim 2] C aoai C ao g 2 Cl ia2 

CL i—l 0,2 — 1 

= - lim [C(0) 4 + (n - l)C(0) 2 C*(oo) 2 + 2(n - 1)C(0) C(oo) 3 + (n - l)(n - 2)C(cx)) 4 l 

2 n— >0 

= ^ [C(0) 4 -C(0) 2 C(oo) 2 - 2C(0) C(oo) 3 + 2C(oo) 4 ] (67) 

We note that the sum produces all not-equivalent diagrams shown of Fig. [TU1 Indeed the first term in the second 
line that follows from a^ = ai = ao corresponds to the first diagram shown in Fig. [TU1 The weight 1 reflects the fact 
that there is only one possible choice ai = a± = aQ. Similarly the second term obtained for a2 = a± =/= ao corresponds 
to the second diagram in Fig. [TUJ There are n — 1 choices that satisfy the constraint a^ = a\ ^ a$. The third term 
follows from either a\ = ao a2 or a\ ^ a§ = <22 and indeed corresponds to the third and fifth diagrams of Fig. 1101 
In both cases there are n — 1 choices that satisfy the constraint. Finally the last term is obtained by taking all three 
indexes ao, a\ and a2 different from each others, and hence [n — V)(n — 2) possible choices. This term corresponds to 
the fourth (or the equivalent sixth) diagram of Fig. [TU1 



B. The self-energy diagrams 

To finalize the calculation of the self-energy diagrams (TT5]) in the limit t — t' — > oo we also need the contribution 
from the empty vertex, the quantity A^(t') in eq. (|32[) . The empty vertex is given by the equilibrium diagrams of 
the associated static theory described by the canonical distribution (TJJ, as a consequence the empty vertex is made 
by the base diagrams used to construct the full- vertex with all lines equal to the equal-time correlation C(0).[37| 

The contribution from any one of such diagram B to the empty vertex can be readily written in the replica formalism: 

JV i-l 

Af(t') = M(B)l[l[K ij C% (68) 

i=l j=0 

where we used the same vertex numbering convention as the full vertex diagrams. The limit n — > is not necessary 
since we are only using the diagonal part of C a b- 
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By using the representations (j66j) and (j68)) the contribution A^ B< _(i, t') of the self-energy diagram with the 1 + r 
full vertex generated by the base diagram Bf and the 1 + r empty vertex generated by the base diagram B e reads in 
the limit t — t' — » oo, see eq. 



Ag Ba (*,*') = M(B f )M(B e )M(B f ,B e 



lim 



l,n 



iV 



/ 4-1 



n n 4 ^ c7 (°°) r n n *$> 



(69) 



.,ajv 2=1 j=0 



i=l j=0 



where M(Bf, B e ) is the symmetry factor generated by the connections between the two vertices. 

The term C(oc) r follows from the lines connecting the 1 + Ef external vertices of Bf to the 1 + E e external vertices 
of B e . Thus by assuming that the replica indexes a, of Bf never take a value equal to that of the replica index b of 
B„ this term can be written as 



C(oo) 



nn^ 

i=0 j=0 



c; 



(70) 



where the rectangular symmetric matrices Jij and rriij give the topology of the connections between Bf and B e , i.e 



1 if the two external vertices are connected or otherwise, while m, 



1,2, 



r) gives the number 



of lines between the two vertices. 

The values of the replica indexes of the two vertices can be made different either assuming b = n + 1 or forcing 
the two sets of replicas to assume different values by adding the projectors (1 — 5 ai b)- This second method is more 
appealing since it maintains the symmetry of replica index, the little price to pay is that now the replica index of 
diagram Bf can take only n — 1 values, so the limit must be changed from n — * to n — > 1. 

Collection all terms we end up with 



A^L (*,*') = M(B f )M(B e 



M (Bf,B e ) lim 



l,n 

E 



N 



S i-1 



nn 



i=lj=0 



3=0 



N e i-1 
i=lj=0 



n ^bb 



(71) 



In Sect. Ill Bl we have seen that the equilibrium dynamic diagrams for T,^(t,t') can be obtained by dividing 
a suitable base diagram with the same topology of the equilibrium dynamic diagrams into two sub-diagrams, one 

leading to Bf and the other to B e . This is precisely the role of the projector n5f (l — ^o«i>) m ec l- |7T|) - The total 
contribution from the self-energy base diagram B is now obtained by considering all possible divisions of B into Bf 
and B e and summing up the result. It is not difficult to realize that all possible distributions of the bare vertices of 
B between Bf and B e can be generated by considering all possible insertions of the projector n£o(l — ^a«&)- ^ s a 

( B) 

consequence the total contribution (t, t') to the self-energy from all equilibrium diagrams generated by the base 

oo reads 



diagram B in the limit t — t' 



X ( B J(t,t>)=M(B) lim 



E 



N+l i-1 

nn^ 

i=l j=0 







(72) 



where we have associated the index with the external vertex at t and the index N + 1 to the external vertex at t' 
so that the total number of bare vertices of B is N + 2. 

The formula has a simple meaning: to obtain the contribution from all diagrams generated by the self-energy base 
diagram B just attach to each vertex of B a replica index a and to each line connecting the vertex a and b the 
correlation function C a b [eq. (|64[) ] . Then sum over all replica indexes from 1 to u keeping the replica index of the 
vertices attached to the t and t' external legs fixed and different from each other. At the end take the limit n — > 1. 

To illustrate the procedure we consider the base diagram B shown in Fig. [3] Equation (f72|) . with the symmetry 
factor M(B) = 1/4, gives 



1 n 

= 4 I™ E C loa, Ca Q a 2 (1 - <5 aoQ2 ) 



- lim [2 C*(0) 2 C(oo) 3 + (n - 2) C(oo) 5 ] 
i[2C(0) 2 C(oo) 3 -C(oo) 5 ]. 



(73) 
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FIG. 11: Fourth order crossed self-energy base diagram. 

The base diagram of Fig. [3] produces the equilibrium dynamic diagrams shown in Fig. [5] whose value is given by eq. 
([20|) , Evaluating the integral for t — t' — > oo we indeed have 

4^*'*') = \f ti d^C{t,s)C{t, s )C{^f + \c{i)fC{^f 
= \ [C(0) 2 - CM 2 ] C(oo) 3 + ic(0) 2 CM 3 

= i[2C(0) 2 C(oo) 3 -C(oo) 5 ]. (74) 

As further example consider the base diagram of Fig. [TT] This diagram has a symmetry factor of 1/2, thus the 
replica calculation gives 



n n 



^ip0 (M ) — 2 j^iJi ^ y ^ y ^aoiil Ca a 2 Caia 2 ^ao<i3 Cai<J3 ^0203(1 ^0003) 



ai— 1 a 2 — 1 

= - lim [2C(0) 4 (7(oo) 3 + nC'(0) 2 C'(oo) 5 +4(n-2)C'(0)C , (oo) 6 + (n-2)(n-3)C'(oo) 7 l 
= i [2C(0) 4 C(oo) 3 + C(0) 2 C(oo) 5 -4C(0)C(oo) 6 + 2C , (oo) 7 ] (75) 
The direct calculation of the equilibrium dynamic diagrams generated by this base diagram is shown in Appendix [B] 

C. The replicated system and the bifurcation equation 

The self-energy base diagram B are the static self-energy diagrams generated by the Hamiltonian H of the associated 
equilibrium static theory described by the canonical distribution (|4|) . Thus from the results discussed so far it follows 
that 

lim Y,00(t,t') = lim E ah , a^b (76) 

t— t'^oa n— >X 

where E a & is the static self-energy of the n times replicated system described by the Hamiltonian Ti n = H[(p a ] 
and the 1RSB correlation function [eq. (f64|) ] 

C ab = <¥W>cq = C(0) 8 ab + C(OO) (1 - 6ab) (77) 

with C(0) = C(t,t) the equal time value of the equilibrium two time correlation function C(t,t') and C(oo) its 
t — t' — > oo limit. The average (• ■ -^q is taken with the canonical probability distribution function 

P " q cxexp(-/3ft„). (78) 

It is now straightforward to derive the bifurcation equation. The correlation C a b and the self-energy E t> are related 
by the Dyson equation 

]T [{CXj - E ac ] C cb = S ab (79) 

c 

where [C ]^ 1 = H"(0)5 ab . The self-energy Y> ab , given by the sum of all 1PI diagrams generated by the vertices of 7i n 
and correlation C ab , has the 1RSB structure: 

E ob = Ei 5 ab + S (1 - 6 ab ). (80) 
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where the diagonal part Ei = E aa is equal to the static self-energy of a single replica, i.e., to the static self-energy 
E^(0) of the original system. By inserting eqs. (|TT|) and ([50]) into the Dyson equation (|T9")) . and summing over the 
replica index c with a 5^ b we have 

[H"(0) - Ei] C(oo) - E C(0) - (n - 2) E C(oo) = (81) 

which, using (|76p . in the n — > 1 limit reduces to 

[H"(0) - S^(0)] C(oo) + [C(oo) - C(0)] E^(oo) = 0. (82) 

Setting now a = b in the Dyson equation (|79|) and taking the limit n — > 1 we end up with 

[ff"(0) - E^(0)] C(0) - (n - l)EoC(co) n = 1 [ff"(0) - E^(0)] C(0) = 1. (83) 
that with eq. (|82p leads back to the bifurcation equation (TH| derived from dynamics. 

V. CONCLUSIONS 

In this paper motivated by the replica formalism developed for the analysis of glassy system, with or without 
quenched disorder, we have shown that the long time limit C(oo) = \im t -t'^cx> C(t,t') of the equilibrium two time 
correlation function C(t, t') can be computed from the n — ► 1 limit of the static of a n-times replicated system with 
a 1RSB structure. In particular we have shown that the Dyson equation of the replicated system leads in the n — > 1 
limit to the bifurcation equation for the ergodicity breaking parameter C(oo) derived from dynamics. The main result 
is the equivalence of the t — t' — > 00 limit of the equilibrium dynamic self-energy H^(t,t') with the n, — > 1 limit of 
the off-diagonal static self-energy E a b of a n-times replicated system with a 1RSB structure. 

The proof is based on the analysis of t — t' — > 00 limit of the self-energy diagrams of the dynamic field theory 
perturbation approach, and follows the following steps. We first discussed the general structure of the equilibrium 
dynamical perturbation diagrammatic expansion of the self-energy E^(t, t 1 ) and shown that this can be written in 
terms of groups of dynamic diagrams obtained by a common base diagram. Next we derived an explicit expression 
valid for any t and t 1 for the contribution of each group. In the limit t — t 1 — > 00 all time integrals can be evaluated 
and for each group one obtains an expression that depends only on C(0) and C(oo). We have proved this result by a 
sum-rule approach and then by a graphical diagrammatic integration method. Next we have shown that the diagrams 
that result from the diagrammatic integration method can be reproduced by introducing "replicas". Moreover by 
choosing for the replicated system a 1RSB scheme we are able to reproduce the correct value of each diagram by 
taking limit n — ► 1 for the number n of replicas. This is proved by showing that in this limit the replica approach 
leads to the same sum-rules derived from dynamics. 

In this paper we focused on the general proof of the equivalence between the replica and the dynamical approaches. 
For this reason we did not discussed any physical consequences of this equivalence. These will be considered in a 
forthcoming paper [21| . 

Finally we observe that the results reported in this paper follow from our studies on the Mode Coupling Theory 
approach to the glass transition. This is why we focused on the self-energy E<^(£, t'). However the method developed 
here can be extended to deal with the t — t' — » 00 limit of other quantities. 
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APPENDIX A: THE LIMIT n -> 

In this Appendix we show that the sum over all replica indices ai with i = 1, . . . , N of any 1 + r base diagram B 
reproduces in the limit n — > the sum rules discussed in Sec. MI Bl 
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Let us denote by 



F B 2) {s h ,s k } = 



Fb{sq; s»!, ■ ■ ■ , Si N _ 2 , Sft, s fc } 

(il,...,ijv-2)GP[(l,...,Af)-(/i : fe)] 



(Al) 



the contribution to As(so) [eq- l|48p ] from all time ordering of s±, . . . , sn with s k the smallest time and s k the next 
smallest one, i.e., s; > Sh > Sfc for any i ^ h,k. Then Ab(sq) can be written as, see eqs. ([50|) and (152j) . 



Ab(so) = J! [4 a) K.**} + i; B ) {**.*fc} 



(/l,fc) 



(A2) 



(h,h) 



We may select in a similar way from the replica sum in eq. (|66p all terms in which the replica indices ah and a/c have 
a value smaller or equal to all others replica indices a^, i ^ h,k. Then by denoting their sum by: 



the equation (|66[) can be written as: 



1,N n 

n e 

i^h.k a,i = l 



N i-1 

n n ^ ' a,i - k - a "^' i ' i: 

i=ij=o 



(A3) 



(A4) 



(h,k) anak 



where the first sum runs over all distinct pairs h, k while the second over ah and a k . The prime ' over the sum sign 
means that only values of a k j t < a?: with i ^ h,k are included into the sum. Finally by decomposing the sums over 
ah and a k as 



E' = E'+ E' + E' 

Qfeflfc a h =a k a h >a k a k >a h 



(A5) 



we end up with 



A B (t) 



hm > 

n-vO ^ 

(/t,fc) 



/ y a.ha k 



(A6) 



It is easy to realize that the three terms correspond to the three terms in the second line of eq. (|A2[) , in the same 
order. The first one is immediate. In the second term a k is always smaller than ah and hence smaller than all ai with 
i ^ k. This is exactly the same structure of the second term of the sum rule (|A2[) . Similarly the third one corresponds 
to the third and last term in (|A2[) . Only the signs are different. The correct signs, and hence the correct second order 
sum rule, are recovered in the limit n — > since in this limit due to the restriction in the sums the second and third 
terms acquire a negative sign from the lines connecting the vertices h and k. 

In a similar way one recovers the third order sum rule. Indeed with a straightforward extension of the notation, 
the third order sum rule, eqs. (|53p and (|57p . can be written as 



A B (so) 



E 



E 



F B 3 '{s ;si : s h ,s k } 
(Ma) P(i,ft,fe) 

|^b 3) { s o; Sl\sh 



(l,h,k) 



si,s k = Sh} 



- F^isilsh = t', s k = s h } - F B *>{s h \s k = t',si = s k } - F B 6 '{s k \si = t', s h = s;} 



.(3) 



(3), 



-F B 3) {si\s h =s u s k = t'} + F ( B 3> {si\s h =t',s k =t'} + F B * i {s h \s l =t',s k =t'} 



r(3) 



(3), 



-F B 3} {s h \s k = s h ,si = t'} + F B 3) {s h \s k =t',si = t'} + F B 6) {s k \s h =t',si =t'} 
- F B 3) {s k \ Sl = s k ,s h = t'} + F B 3) {s k \ Sl = t', s h = t'} + F B 3) { Sl \s k = t', s h = t'} 



.(3) 



(3), 



(A7) 
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where 

F^' {si, s h , s k } = FB{so- ) s il ,...,Si N _ 3) 8i,ah,Sk} (A8) 

(i 1 ,...,is-3)eP[(l,...,N)-(l,h,k)] 

On the other hand by selecting three replica indexes, the replica expression (|66p leads to 

M*)=l™ E E^fU ( A9 ) 

(l,h,k) o.ia h a k 

where 

II II ^ G ^., , aj.ft.fc < Oi, i ^ Z, ft, ft (A10) 

i=l j=0 

Now by decomposing the sum over the three replica indexes as 

£' = £' 

ajajjOfe a l =a h =a k 

+ E'+ E'+ E' 

rai>Oh=ajs a fc >a fc =a; a fe >ai=a t 

+ E'+ E'+ E' 

ai=ah>afc ai>a h >a k a h >a t >a k 

+ E'+ E'+ E' 

a h =a k >ai a h >a k >ai a k >a h >ai 

+ E'+ E'+ E' (aid 

ifc=ai>ah a k >ai>a h a t >a k >a h 

one identifies all the terms of the third order sum rule (|A7|) . The third order sum rule (|A7|) is recovered in the 
limit n — > since each inequality in the sum gives in this limit to a "minus" sign which, when combined together, 
reproduces the correct sings. 

From these two examples it should be clear that the correspondence between the sum rules derived with the two 
formulations, dynamic with integrals and replica with sum over integers, can be extended to sum rules of any order. 
Indeed the multiple sums in eqs. (|48[) and (|66[) can be always decomposed in the same way by peaking up the same 
group of vertices. Moreover the lower limit of the integrals in the dynamical formulation is always associated with an 
inequality in the correspondent sum of the replica formulation. Then, since in the limit n — > each inequality gives a 
"minus" sign, it is clear that in the n — ► limit both formulations lead to the same sum rules and hence to the same 
expression for As(i), once the two point correlation function C a b of the replicated system is identified with (|64[) . 

APPENDIX B: THE FOURTH ORDER CROSSED DIAGRAMS 

In this Appendix we sketch the calculation of the equilibrium dynamic diagrams generated by the fourth order 
self-energy base diagram of Fig. [11] in the limit t — t' — > oo. The equilibrium dynamic diagrams generated by the base 
diagram of Fig. [TTJare shown in Fig. [TJ1 see Section Hi Al and lead to: 

Sj? (*>*') = j t ds j t rfs ' G ( t I s ) c '(^ s, )G(s,s , )C(s,s , )C(s,<')C(s',i , )C(^i') 
+ J ds' J dsC(t,s)G(t,s')G(s',s)C(s',s)C(s,t')C(s',t')C(t,t') 

+ \J ds J ds ' G (^ s ) G (^ s ')C(s,s') 2 C(s^')C(s',<')C(^^) 

1 /"* 

+ - / dsG(t,s)C(s,t') 3 C(t,t') 2 C(t',t') 

2 Jt' 



F(3) 

aio.h,a k 



M(B) 



1,N 



n e 
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FIG. 12: Fourth order crossed equilibrium dynamic diagrams. 



+1 f ds'G(t,s')C(s',t'fC(t,t/) 2 C(t',t') 
2 Jt' 

+ l -c{t,t'fc{t'X) A - 

The factors 1/2 are the symmetry factor of the diagram. 
The third integral can be split in two by using the identity 



(Bl) 



ft ft ft f-S ft ps' 

/ ds ds' = ds ds' + / ds' / ds (B2) 
Jt' Jt' Jt' Jt' Jt' Jt' 

By using now the FDT relation G(s, s r ) — d s 'C(s, s') the first and half of the third integral can be written as: 

\ I ds [ ds'd s C(t,s)C(t,s')d s >C(s,s') 2 C(s,t')C(s',t')C(t,t') 
2 Jt' Jt' 

+ \j t ds / ds ' d s c ( t ^) d s'C{t^')C{s,s'fC{s,t')C{s\t')C{t 1 t') 



2Jt 



ds / ds' d s C(t,s)d s/ [C(t,s')C(s,s') 2 ]C(s,t')C(s' ,t')C(t,t') 



(B3) 



In the limit t — t' — > oo we can replace each of the last three correlation functions by C(oo) since the integrand is 
different from zero only if s ~ s' ~ t ^> t' . The integral over s 1 can now be easily done and one gets 



ic(0) 2 C(oo) 3 f dsd s C{t,s)C{t,s)-\c{oof ( dsd s C(t,s) 
2 Jt' 2 J t i 

= ic(0) 2 C(oo) 3 [C(0) 2 - CM 2 ] \c{<x>f [C(0) - CM] 



(B4) 
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The other half of the third integral can be combined with the second integral, the one from diagram (6), and a 
similar manipulation leads again to the result (|B4|) . This is not unexpected since diagram (a) and (6) can be changed 
one into the other by exchanging s and s' . 

The fourth and fifth integrals, i.e., those from diagrams (d) and (e), can be readily done with the help of the FDT 
relation and in both cases one obtains for t — t' — > oo: 

ic(0)C(oo) 5 [C(0) - C(oo)] . (B5) 

Collecting all contributions we ends up with 

lim 4 4 - c) (i, t') = l\2 C(0) 4 C(oo) 3 + C(0) 2 C(oo) 5 - 4 C(0) C(oo) 6 + 2 C{oo) 7 } (B6) 

as found from the replica calculation, see eq. ([75)) . 



APPENDIX C: DERIVATION OF EQ. © 



Our aim is to study the correlations associated with the stochastic process (TTJ). Instead of working directly with 
eq. (fTJ) it is more convenient to construct a generating functional from which correlations can be obtained. Following 
the standard procedure of field theory [29| one introduces an external time-dependent source J(t) and defines the 
generating functional [H[ 



/• r roo r /-ex 

Z\J\=M I VtpUrjV[<po] 5(ip — ifn) exp — / dtJtp exp — / 

J I J t J L J t 



(CI) 



where (p v is the solution of stochastic the eq. ([T]) for a given realization of the stochastic field r/(t) and initial condition 
(po = tp(to). The initial condition is assigned with the probability "Pk^oli that we take equal to the equilibrium 
probability distribution ([4]). Finally J\f is a normalizing constant. The S- function stands for 



6(<fi - tp v ) = S 



at 5tp(t) 



v 



det 


Sr] 




Sip 



(C2) 



where the factor det \Sr]/Sip\ is the Jacobian of the transformation rj — > <p, that in the Ito calculus is equal to one. 

By using the integral representation of the delta function with the help of the auxiliary hat-field tp, and performing 
the integral over the stochastic field 77, a straightforward algebra leads to 



Z[J] = J\f J e^'^-'C dtJip (C3) 

where S'[ip, (p\ is given by eq. (|18[) and includes the contribution from the initial equilibrium distribution ([4|). By 
adding a second external time-dependent source J(t) coupled to (p all correlations and responses can be obtained from 
differentiation of the generating functional 



Z[J,J] =Af I VipVip e 
Define now the 2x2 correlation matrix 

(t,t>) g v e(t,t?) 



-S'[tp,<p\- f°° dt [J v +Jip] 



G(t,t') 



(t,t r ) g^(t,f) 



(mvV)) mm 1 )) 



C{t,t') G{t,t') 
G(t',t) C(t,t') 



(04) 



(05) 



The average is taken with the dynamical functional S f [<p,<p\. We assumed for convenience (ip) = (<p) = 0. The 
correlation matrix Q(t,t') is solution of the Dyson equation 



/>oo 

/ ds V[G-i(«,s)-£a 7 (M)] Q 1 p(s,t') = 5 a p5(t-t') 

J t„ 



(06) 
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where the Greek indexes run over the two values (ip,0), and the matrix G a p(t,f) is the "free" correlation matrix 
obtained from the second order term of the expansion of S'[ip, (p\ about (ip) = (ip) = 0. A simple calculation gives 

/ H"[0] 6{t - to) 6(t - t') {-d t + H"[0]) S(t - t') 
G-\t,t)=\ | (C7) 

\ (d t +H"[0\) 8{t-t') -2S(t-t') 

where H"[0] = 5 2 H[ip]/Sip 2 \ =Q . 

The equation for the correlation function C(t, t') follows now from the Dyson equation (|C6[) by choosing a = <p and 

= <p: 

[d t + H"[0]] C(t, t') - 2G(t', t)- [ ds E^(t, s) C(s, t') - f ds E^(t, s) G(t' , s) = (C8) 



In equilibrium FDT holds, then inserting the FDT relations 

G(t,t') = -d(t-t')d t c(t,t'), E 0v (t,t') = -e(t-t')d t ^(t,t') (C9) 

into eq. (1C8|) and integrating by parts a straightforward algebra leads to eq. ((9]) of the main text. 
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